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We demonstrate optimal state estimation for a cavity optomechanical system through Kalman 
filtering. By taking into account nontrivial experimental noise sources, such as colored laser noise 
and spurious mechanical modes, we implement a realistic state-space model. This allows us to 
obtain the conditional system state, i.e., conditioned on previous measurements, with minimal 
least-square estimation error. We apply this method for estimating the mechanical state, as well 
as optomechanical correlations both in the weak and strong coupling regime. The application of 
the Kalman filter is an important next step for achieving real-time optimal (classical and quantum) 
control of cavity optomechanical systems. 


Introduction .—State estimation is a crucial task at the 
heart of control theory, both in the classical [1] and in the 
quantum domain [2]. For Gaussian systems, real-time 
state estimation can be achieved in an optimal manner 
using Kalman-Bucy filtering [3, 4]. Since many physi¬ 
cal systems are approximately Gaussian, Kalman filter¬ 
ing has been successfully implemented for a broad range 
of uses, for example for navigation and tracking in aero¬ 
nautics (including the Apollo project and the Global Po¬ 
sitioning System GPS) [5], as well as in the physical sci¬ 
ences, such as for suspension noise cancellation in gravita¬ 
tional wave detection [6] , Heisenberg limited atomic mag- 
netometry [7] or quantum-enhanced optical-phase track¬ 
ing [8-11]. In this Letter we introduce a new domain of 
applications by implementing Kalman filtering for cav¬ 
ity optomechanical systems. These systems represent 
a versatile light-matter interface in which optomechani¬ 
cal interactions inside optical or microwave cavities allow 
control over optical and mechanical degrees of freedom. 
While the first investigations go back to the late 1960s 
in the context of gravitational wave detectors [12, 13], 
it is only the last few years that have seen the devel¬ 
opment of a completely new generation of micro- and 
nano-optomechanical solid-state devices with fast grow¬ 
ing application areas from classical sensing to quantum 
information processing [14]. 

State estimation of a cavity-optomechanical system in 
real time is key for optimal state control and verification. 
The outstanding challenge is to obtain reliable informa¬ 
tion on the mechanical subsystem. In optomechanics, 
this is done through an optical cavity field, which im¬ 
poses both additional noise and dynamical back-action 
effects that have to be taken into account. Until now. 
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reconstructions of the mechanical dynamics have focused 
either on statistical properties [15, 16] or, for real-time 
reconstructions, on regimes of sufficiently weak coupling 
and negligible dynamical back-action effects [17-20]. The 
information obtained from real-time estimation about 
the mechanical quadratures can be used for active feed¬ 
back control of the mechanical resonator [19-21]. How¬ 
ever, the validity of these reconstruction schemes breaks 
down when either coupling strength, dynamical back- 
action effects or noise become strong. Our Kalman fil¬ 
tering approach overcomes this limitation and allows us 
to demonstrate real-time optimal state estimation for 
cavity optomechanical systems operating in arbitrary 
parameter regimes. From a quantum physics perspec¬ 
tive, the Kalman filter solves the stochastic Schrodinger 
equation—a stochastic, nonlinear generalization of the 
Schrodinger equation—which is the canonical way to de¬ 
scribe quantum systems subject to a continuous mea¬ 
surement via coupling to electromagnetic fields [2, 22- 
24]. These concepts and their application to mechanical 
systems have been the subject of extensive theoretical 
research [25-29], but no experiments in the context of 
cavity-optomechanics have been conducted so far. 

Kalman Filter .—In quantum theory, just as in classi¬ 
cal theories, a continuously observed system can be de¬ 
scribed by a conditional state [2], i. e., a state that incor¬ 
porates the total amount of knowledge that an observer 
has extracted from her set of measurements. Discarding 
this knowledge yields the unconditional state, which is 
an incoherent mixture of all possible conditional states. 
Our goal is to find the (multipartite) conditional states 
of the full cavity-optomechanical system including me¬ 
chanical and optical subsystems. We restrict ourselves to 
Gaussian dynamics and measurements, which is a valid 
assumption for the existing realizations of optomechan¬ 
ical systems [14]. For this case, it has been shown [30] 
that the problem of finding the conditional state can be 
mapped to a classical estimation problem, which is solved 
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by a Kalman filter. It produces a real-time state esti¬ 
mate from a continuous measurement trajectory, which 
is optimal in the sense of minimizing the mean square es¬ 
timation error. We describe the system by the following 
(linear) state-space model 

= + (la) 

Zt = CtXt+Vt, (lb) 

where x^ is a state vector in some appropriately chosen 
state space (e. g., the phase space of a harmonic oscilla¬ 
tor), Zt is the outcome of a linear measurement on the 
system, and and describe process and measure¬ 
ment noise, respectively. Both and are assumed to 
be zero-mean Gaussian white-noise processes, which obey 
Re(E[wtwJ]) = W6{t — s) and Re(E[vtvJ]) = V6{t — s), 
where 6 is the Dirac (i-function, E['] denotes the expec¬ 
tation value with respect to the initial probability dis¬ 
tribution describing system and noise, and Re(-) the 
real part Process and measurement noise may be 
correlated, which is described by the cross-correlations 
Re(E[wtvJ]) = MS{t — s). At and Ct are real, matrix¬ 
valued functions that parametrize the state-space model 
and are fixed by the physical model of the system and 
measurement process. We aim to find the estimate x^ 
of Xt that minimizes the mean square estimation error 
E[||xt — Xtip] at a time t when taking into account the 
measurement results {z^ : 0 < s < t}. This yields 
Xt = E[xt|{z 5 : 0 < s < t}], i.e., the conditional ex¬ 
pectation value of Xt given the full measurement record. 
Evaluating this for system (1) yields the time-continuous 
version of the Kalman filter [4] 

ict = At±t + Kt{zt-Ct±t), (2a) 

Pt = AtPt + PtAj + W- KtVKj, (2b) 

where = {PtCj+M)V-'^ is the so-called Kalman gain 
and Pt = Re(E[(x( - X()(x( - x*)'’’]) is the estimation- 
error covariance. 

We can ascribe a quantum theoretical meaning to x^ 
and Pt by associating x with the Schrodinger operators 
that describe the quantum system, x^ with the corre¬ 
sponding Heisenberg operators that evolve under (la) 
(their Heisenberg-Langevin equations), and z^ with an 
operator-valued output process [24, 30]. If pt is the Gaus¬ 
sian state conditioned on a continuous measurement of 
zt, we have Xt = tr{xpt} and Pt = Re(tr {xx^pt}) — 
x^x^, i. e., the symmetrized covariance matrix of x^ with 
respect to p. In other words, the conditional Gaussian 


^ Taking the real part of the covariance matrices that describe the 
noise processes is only necessary for quantum processes due to 
their non-commutative nature. Although it is not necessary for 
classical systems, we choose this explicitly real form for the sake 
of a consistent presentation. 



FIG. 1. (color online) Kalman filter for cavity-optomechanical 
systems, (a) Working principle of the Kalman filter: (i) The 
conditional state is depicted by the green phase-space ellipse, 
which (ii) evolves in time according to the system dynamics, 
(iii) After a time dt a Bayesian update is applied based on 
the measurement outcome to find the new conditional state. 
This procedure minimizes the mean-square estimation error, 
which makes the Kalman filter optimal for real-time state 
estimation, (b) Schematic of the experiment: The optome¬ 
chanical cavity is driven by two laser beams each of which 
carry amplitude- and phase-noise. The mechanical motion is 
typically driven by Brownian noise. After their interaction 
with the cavity the optical fields are detected by two inde¬ 
pendent homodyne measurements (signals Zd and Zr), which 
themselves are subject to optical losses and noise. Building 
an accurate Kalman filter requires appropriate modeling of 
all relevant noise sources. 

state pt is parametrized by x^ and Pt ^ • By averaging 
over all possible trajectories of x^ we recover the un¬ 
conditional state, whose covariance matrix Re(E[xtX^]) 
we can extract from the estimated data by noting that 
E[(xt —xt)x^] = 0, and thus Re(E[xtX^]) = Pt+lE[xtX^] 
[2, 31], 

The Kalman filter equations (2) describe how the con¬ 
ditional state is iteratively updated (Fig. la). First, the 
estimate Xt and the covariance Pt are propagated for an 
infinitesimal time interval dt [first term in (2a) and first 
three terms in (2b)] according to the state-space model 
(la). Second, the measurement outcome is incorporated 
as a Bayesian update that corrects the value of the esti¬ 
mate Xt and contracts the covariance ellipse [last terms in 
(2a) and (2b)]. The updated values are again propagated 
by dt and the procedure is repeated. 

The model .—We consider a typical cavity- 

optomechanical architecture (Fig. lb), in which a 


^ One can also adopt a quantum-optical interpretation of xt. For¬ 
mally integrating (and assuming vanishing initial conditions) 
gives xt = f—oo s)zsds with an integral kernel K depend¬ 
ing on At, Ct and Kf. Thus xt is formally equivalent to an 
(unnormalized) bosonic mode extracted from the output process 
zt. In a quantum-optical setting this could be for example a 
temporal light mode extracted from the output light of a cavity. 
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Fabry-Perot cavity (resonance frequency Uc), coupled 
to a single mechanical mode is driven by two laser 
fields (at frequencies ^^o,r)- The “resonant” beam 

{ujr = 00c) acts as a weak probe of the cavity length 
to stabilize the laser frequency with respect to the 
cavity resonance; the “detuned” beam {ood ^ ooc) induces 
dynamical back-action effects, e. g., for laser cooling. 
This captures all relevant scenarios applied in typical 
optomechanics experiments. The mechanical element 
has a resonance frequency oOm and energy damping rate 
jrn- Both cavity modes exhibit decay at a (half width 
at half maximum) rate k = ki K 2 , where ki describes 
the input coupler and tv 2 accounts for spurious photon 
losses. The system is described by the (linearized) 
quantum Langevin equations [32-35] 

q = uJmP (3a) 

P = - ImP + XI 5 i{cos OiXi - sin OiVi) + (3b) 

i=r,d 

Xi = -KXi + Ajj/i + Qi smBiq + 

+ + 2^/K^S/3i + |Q!o,i| sinOi^i, (3c) 

m = -K-Vi - ^iXi + gi cos 0iq + 

+ + |ao,i| cos6>j0j, (3d) 

where p {[q^p] = i) describe position and momen¬ 
tum of the mirror, and pk with [xi^pk] = for 
l^k G {r, d} respectively denote the amplitude and the 
phase quadrature of the cavity modes of the resonant 
and detuned beam. The optomechanical coupling to 
cavity mode i is given by pi = \/2^o|<ao,i| with = 
!hjjjQ^i!{n -1- iA^), where qq is the single-photon 
coupling strength. Pi is the corresponding driving laser 
power, and = co’o,i — oOc is the detuning of the re¬ 
spective driving laser (at ooo^i) with respect to the cavity 
resonance frequency (cJc). The coupling of the mechan¬ 
ics to a thermal bath is modeled by a self-adjoint noise 
term ^ with {^{t)^{s) + ^(s)^(t)) = 27 ^( 2 n + l)6{t - s) 
and n ~ kBT/hoJm (the mean occupation number of the 
bath at temperature T). Optical shot noise is denoted by 
VT with variances {xf^{t)x^j^{s)) = {yl^{t)y^^{s)) = 
^6ij5{t — s). Terms proportional to 6Pi and pi describe 
classical amplitude and phase noise of the driving lasers 
[33-35]. 

Homodyne detection is used to independently measure 
the generalized quadratures Zi of the reflected optical 
modes (Fig. lb). The cavity input-output relations yield 

z'i = Xi -h x^P + 5Pi) cos pi 

+ (V2^2/i + 2 /‘") sm ipi, (4) 

where 5Pi describes classical amplitude noise. We model 
optical losses and inefficient detection as beam-splitter 


^ The generalization to several mechanical modes is straight¬ 
forward, and, in fact, is included in the full model of our system. 


losses parametrized by y. The measured quantities are 
rescaled to Zi = ^/1 — y z^ + ^Jy zf ^, where zf^ describes 
additional quantum noise independent of x^p and yp^ i. e., 
{xp{t)zp{s)) = {yp(t)zp{s)) = 0. Defining the vec- 
tors xt = {q{t),p{t),Xd{t),yd{t),Xr{t),yr{t))'^ and zj = 
{zd{t)^ Zr{t)) , equations (3) and (4) can be rewritten in 
the compact form (1). 

Contrary to the idealizing assumptions made above, 
many of the noise sources in an actual experiment 
are frequency-dependent, here the laser amplitude- and 
phase-noise. This needs to be taken into account by 
properly extending the state space model. We incorpo¬ 
rate three types of laser noise: (i) broadband laser noise 
originating from the laser itself, (ii) narrow-band Pound- 
Drever-Hall phase modulation in the resonant beam re¬ 
quired for locking the laser to the cavity frequency, and 
(iii) narrow-band laser noise originating from the feed¬ 
back loop of the laser lock. Each of these noise sources 
is experimentally characterized and is modeled indepen¬ 
dently to match the overall spectral characteristics (see 
Appendix C). Furthermore, we extend the state space 
model to incorporate higher-order mechanical modes. 

Measurements and innovations. —We use the recorded 
homodyne signals Zi as input to the Kalman filter for esti¬ 
mation of the optomechanical state, which is done offline. 
Fig. 2 shows a 2/is trace of the detector signals (corre¬ 
sponding to 100 sample points), along with the optimal 
measurement prediction. The prediction shows excellent 
qualitative agreement with the measured data both in 
the weak {gd < n) and in the strong coupling regime 
idd > ^)- Quantitatively, the validity of the estimation 
is assessed by the innovation sequence Uf = Zf — Ct^t^ 
i. e., the difference between the predicted measurement 
Zf = Cf'^t of tfio Kalman filter and the actual measure¬ 
ment outcome z^. For an optimally working filter, 
must be a Gaussian zero-mean white noise process with 
a variance given by ’Eliyfiyp] = CtPfCp + V. We use this 
fact to fine-tune model parameters starting from their 
independently determined values. The statistics of Uf of 
the resulting Kalman filter closely matches these criteria, 
hence demonstrating the accuracy of the filter (Fig. 2; see 
also Appendix D for further statistical analysis). 

Estimation of optomechanieal quadratures. —Kalman 
filtering provides direct, real-time access both to the opti¬ 
cal intracavity quadratures and to the mechanical degree 
of freedom in a cavity-optomechanical system (Fig. 3a). 
In the weak coupling regime, the thermally driven me¬ 
chanical motion and its coupling to the optical intra¬ 
cavity fields is visible. Clearly, the mechanical motion 
modulates both quadratures Xd^ Pd of tho detuned beam 
[14], which couples to the mirror via the optomechani¬ 
cal beam-splitter interaction. The situation is different 
for the resonant beam, whose amplitude quadrature Xr 
contains shot noise only, while its phase quadrature Pr 
couples to the mechanical position. 

The phase space representation captures the essence of 
Kalman filtering. The estimated mechanical quadratures 
rotate in phase space (Fig. 3b). Their probability distri- 
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weak coupling (gd=0.2 k) 


strong coupling (gd=1.68 k) 



Time (ps) 



Time (ps) 


FIG. 2. (color online) Measurement signals and Kalman- 
filter predictions. Shown are the measurement signals of 
the two homodyne detections of the detuned and resonant 
beam, Zd{t) (red dots) and Zr{t) (blue dots), respectively, and 
their Kalman-filter predictions Zi(t) (gray line) both for the 
weak (left) and the strong coupling regime (right). Error 
bars of the prediction (±2(7) are indicated by the width of 
the gray line. Kalman filter innovations z^z(t) = Zi(t) — Zi(t) 
are plotted below each data set and demonstrate the accu¬ 
racy of the implemented filter. To assess the performance 
of the filter, we calculate the fraction of normalized innova¬ 
tions that are contained in a two-sided 95% confidence region 
(±2(7, indicated by the lines) of a zero-mean Gaussian distri¬ 
bution (besides each plot). The experiment was performed 
at room temperature with a micromechanical oscillator of 
ujm = 27r X 1.278 MHz, 7 ^ = 27r x 265 Hz, and optomechan¬ 
ical parameters k = 0.34 ct;m, go = 27 t x 7.7 Hz (for details 
see Appendix B). For both coupling strengths of the detuned 
beam (A^ = curn) we use (fd ~ 0, = 0, pr = 0.2and 

(fr = 7t/2. Note that the fast oscillation of Zr{t), which is 
due to the 20 MHz Pound-Drever-Hall phase modulation for 
frequency locking, is taken into account by the Kalman filter. 


bution along each mechanical quadrature is shown as his¬ 
togram besides each axis and demonstrates the Gaussian 
nature of the micromirror motion. We compare the un¬ 
certainty ellipse of the unconditional (dashed line) and 
conditional (solid line) mechanical state, i. e., the area 
in which we expect with 95% probability to find the me¬ 
chanical quadratures. For a purely thermal state the area 
of the unconditional ellipse is proportional to the ther¬ 
mal occupation number n. In the weak coupling regime, 
the information provided by the measurement update 
leads to a clear reduction in the uncertainty (factor of 
27 in effective temperature), which is the optimal one 
for the given coupling strength In the strong coupling 
regime, laser cooling has already significantly diminished 
the thermally induced uncertainty of the unconditional 


Recall that the optimality of the Kalman filter ensures that the 
conditional state uncertainty is minimized for a given coupling 
strength. 



FIG. 3. (color online) Estimating optomechanical quadra¬ 
tures. (a) Shown are Kalman-filter real-time estimates for 
the optical amplitude (xi, straight line) and phase (y^, dashed 
line) quadrature of the detuned (red, top) and resonant (blue, 
middle) beam along with the mechanical position {q, straight 
green line) and momentum quadrature (p, dashed green line) 
of the optomechanical system for the weak and strong cou¬ 
pling regime. The mechanical (b, green line) and optome¬ 
chanical (c, gray line) phase space trajectories are estimated 
over a period of lOO/xs. A histogram along each quadrature 
is shown as side panel and estimated over 10 ms. The uncer¬ 
tainty ellipse of the unconditional (conditional) state is shown 
as dashed (straight) line. Note that the length of the shown 
trajectory is not sufhcient to adequately represent the state’s 
statistics. All units are given in terms of quadrature zero- 
point fluctuations (zpf). For our experimental parameters 
Qzpi = 2.73 X 10“^®m, pzpf = 3.87 x 10“^^kgm/s. 


state. In addition, the cavity dynamics introduces a no¬ 
table ellipticity in the phase-space distribution [36]. The 
conditional state uncertainty is similar to the weak cou¬ 
pling situation. This is because for technical reasons the 
signal power at the homo dyne detectors was kept con¬ 
stant for both coupling strengths, which means that the 
stronger detuned optical drive beam does not provide 
more information on the system state. 

Figure 3c shows real-time estimates of the optome¬ 
chanical correlations between mechanical position and 
the phase quadrature of the resonant beam. Analogous 
to the mechanical phase space, the conditional state un¬ 
certainties are strongly reduced, reflecting the real-time 
information gain on the optomechanical correlations. 

Conclusion .—We have successfully implemented 
Kalman filtering for optimal state estimation of cavity 
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optomechanical systems. Its accuracy crucially relies on 
an accurate state space model of the specific experiment. 
The applications of this method in the domain of 
optomechanics are manifold. For example, Kalman 
filtering enables mechanical feedback control in the 
quantum regime. While in this work we operate the 
filter offline, its real time application in the frequency 
range investigated here is feasible using current field 
programmable gate array hardware (see Appendix 
E). The optimality of the filter guarantees that the 
reduction in conditional state uncertainty corresponds 
to the maximal cooling one can achieve through ac¬ 
tive feedback at this specific coupling strength. As a 
consequence, ground state cooling is readily achievable 
by combining Kalman filtering with measurements in 
the strong cooperativity regime [31]. This regime has 
been reached in current experiments [21, 37-39]. In our 
case, it requires cryogenic cooling of the mechanical 
environment to 300 mK and quality factors above 10^. 
As another example, mechanical sensing requires precise 
knowledge of the system dynamics in the absence of 
the external impetus, which is equivalent to the task of 
implementing the optimal estimator, i. e., the Kalman 
filter. The same is true for the task of characterizing 
or reconstructing an optomechanical quantum state (for 
example in terms of entanglement), where the relevant 
information is often encoded in the covariance matrix 
Pf. One fascinating prospect there is the generation 
of entanglement of macroscopic test masses through 
measurement [29, 40]. In summary, Kalman filtering 
adds a significant performance advantage for classical 
and quantum control of cavity optomechanical systems. 
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Appendix A: State space models 

1. Systems driven by white noise 

State space models and Kalman filters based on them 
are widely used in classical signal processing [41-43]. 
Here we apply them to an optomechanical experiment in 
a way that remains applicable for quantum experiments. 

Consider a classical Gaussian system with state vec¬ 
tor Xt that is continuously monitored by linear measure¬ 
ments with outcomes z^. In general, both the evolution 
of the state vector as well as the measurement are af¬ 
fected by stochastic noise called process noise and 
measurement noise v^, respectively. Additionally, both 
the state evolution as well as the measurement may be 
subject to a deterministic input u {e.g. a known force). 
Then, the joint (stochastic) evolution of state vector 
and measurement outcomes Zt are described by a state 
space model of the form 

xt = AtXt + BtWt + Lt^t (Ala) 

Zt = CtXt + DtWt + vt. (Alb) 

Here, At^Bt^Ct^Dt and Lt are (potentially time- 
dependent) matrices whose naming is given in Tab. H. To 
completely specify the state space model (Al), we also 
have to specify the statistical properties of the noise. If 
process noise and measurement noise are white and 
Gaussian this amounts to specifying covariance matrices 
W and V for process and measurement noise, respec¬ 
tively, as well as their correlation matrix M. 

2. Colored noise and linear filters 

The Kalman filter relies on a state space model of the 
form (Al) for which the process and measurement noise 
are Gaussian and white. Our experiment, however, is 
subject to colored Gaussian noise, namely classical am¬ 
plitude and phase noise of our laser. But we can extend 
the state space model of the actual optomechanical sys¬ 
tem to an equivalent larger state space model, which is 
only driven by white noise, as described in the following. 
The Kalman filter can then be applied to this extended 
state space model. 

Consider an n-dimensional (colored) noise process ^t 
whose time evolution is described by ^t = + Ct 5 where 

(t is white noise. Its spectrum is given by the rational 
function S^(u;) = where 

H{uj) = —{iuj + F)~^ is the transfer function of the pro¬ 
cess, W is the covariance matrix of ("t, and deg{q^) = n. 
Consider now a system driven by the colored noise pro¬ 
cess he., of the form 

Xt = AfXt + BfUt P Wt + Ct- (A2) 

To model this system by an equivalent state space model 
driven by white noise only, we can define yt = {xj 
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and extend the state-space model (A2) to 



We use such a state space model extension to incorporate 
classical, non-white laser noise. For our experiment, we 
model noise of narrow Lorentzian line-shape as well as 
broadband colored laser noise. 

In a similar way, we can model the effect of an elec¬ 
tronic or digital filter of order n on an input signal 
Uf. In this case, (t is related to the input signal by 
(t = Gut^ and the output is given by + Dut 

(where G, H and D are matrices appropriately cho¬ 
sen to describe the specific filter). Its transfer func¬ 
tion, relating the output to the input signal, is given by 
= yt{s)/ut{s) = H + D with s = iw. 

Appendix B: Experiment 
1. Setup 

The experimental setup is sketched in Fig. 4. Our 
optomechanical cavity is a Fabry-Perot cavity with 
linewidth (half width at half maximum, HWHM) k 
440 kHz and free spectral range (FSR) 15 GHz. Its 
end mirror is a mechanical oscillator with frequency 
27r • 1.278 • 10^ Hz and linewidth (full width at 
half maximum, FWHM) 7^ 27r • 265 Hz. The oscil¬ 

lator consists of a doubly-clamped SiN bridge with a 
Ta205/Si02 Bragg mirror on top. The optomechanical 
cavity is placed inside vacuum {p 5 • 10“^mbar); all 
measurements are taken at room temperature. 

The mechanical oscillator interacts with two opti¬ 
cal cavity modes in orthogonal polarizations and sepa¬ 
rated in frequency by one FSR. These two cavity modes 
are driven by two laser beams, one resonant and one 
red-detuned by approximately the mechanical frequency. 
Both beams are derived from a single NdiYAG laser with 
a wavelength of 1064nm as described in more detail in 
Fig. 4. The laser is locked to the optomechanical cavity 
using a Pound-Drever-Hall (PDH) lock in the resonant 
beam. 

The reflected beams (the signal beams) are split off 
using wave plates and Faraday rotators and directed to 
two separate homodyne detection setups. The local os¬ 
cillators for the homodyne detection are derived from the 
incoming beams. The path lengths of signal beams and 
local oscillators are carefully matched to ensure that sig¬ 
nals and local oscillators always have the same instanta¬ 
neous frequency such that the laser phase noise does not 
affect the homodyne detection. The relative phase of lo¬ 
cal oscillator and reflected (signal) beams can be locked 
using feedback to piezo-driven mirrors (PZT) in the lo¬ 
cal oscillator beam paths. This allows, for each returning 
signal beam, to measure an arbitrary generalized quadra¬ 
ture. For the detuned beam, the phase pd is scanned in 


time between 0 and 7r/2, whereas for the resonant beam, 
the phase is fixed to 7r/2. 

2. Classical laser noise 

a. Laser field 

We describe the extra-cavity laser fields of each beam 
by a displaced coherent state 1/3(t)) with a fluctuating 
coherent amplitude P{t) (note that we suppress the mode 
indices d and r here and in the following). In a frame 
rotating at the laser frequency, we can write the coherent 
amplitude as 

/3(i) = (/3o + W)e-'^W, (Bl) 

where /3(t) is a classical complex random variable with 
expectation value E[/3(t)] = /3o G R. The amplitude 
and phase fluctuations (j){t) are assumed to be real, 

small (E[6/3‘^] /3o,E[(/)^] ^ 1 ), uncorrelated (^[6^(1)] = 

K[613] IE[0]) and zero-mean (E[(i/3] = E[0] =0). Here and 
in the following, we use 

E[/(/?)] = [ p{/3) /(/3) 

Jc 

to denote the classical expectation value of f{f3) and p{/3) 
for the probability density function of (3 in the complex 
plane. Taking the classical fluctuations into account, we 
must therefore describe the extra-cavity laser fields by a 
mixed state 

P = J\P) {P\ ■ 

The expectation value of an operator O in state p is then 
defined as 

{0)=tv{pO} = I d2/?p(/3) (/3|0|/?). 

Assuming that we are only dealing with wide-sense sta¬ 
tionary random processes, the noise power spectrum 
Soo{^) of an operator-valued noise process 0{t) is de¬ 
fined as [44] 

^oo(c^)= [ {0{t)0{0)). 

Jn 

To calculate noise power spectra of photocurrents we 
need to use the following properties of the annihila¬ 
tion and creation operators 6, of the extra-cavity field 
modes 

Kt) \l3{t)) = p{t) |/?(t)) 

(/3(i)|6(t)t = P{t)* W{t)\ 

where S{t) denotes the Dirac delta distribution and the 
operators 5, have units of \/Hz. 
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FIG. 4. (a) Experimental setup. Laser light at 1064 nm is split into two beams. One of the beams (blue) is sent to two free 
space electro-optical modulators (EOMs), which add phase modulations at 1.45 MHz and 20 MHz. The modulation at 1.45 MHz 
is used for calibration of the optomechanical signal. The modulation at 20 MHz is used for locking the laser frequency to the 
optomechanical cavity (OMC) using a Pound-Drever-Hah (PDH) lock. This beam is therefore always resonant with the OMC 
and we refer to it as the “resonant beam”. 

The second beam is sent to a hber EOM which creates strong sidebands at approximately 15 GHz. One of the 1st order 
sidebands is used, whereas the carrier and all other sidebands are blocked by an optical hlter consisting of a volume holographic 
grating and a broadband hlter cavity. Hence, the resulting beam (red) is detuned with respect to the resonant beam (blue) by 
~ 15GHz which is approximately the free spectral range (ESR) of the OMG. This beam is therefore close to another longitudinal 
resonance of the OMG but can be detuned by an additional detuning A by changing the modulation frequency. 

The resonant and the detuned beam are combined on a polarizing beam splitter (PBS) into a single spatial mode but orthogonal 
polarizations. This mode is mode-matched to the OMG, which is located inside a vacuum environment. The returning beams 
(the signals) are separated from the incoming beams using optical circulators (half-wave plates, Earaday rotators and PBS) and 
sent to two homodyne detection setups. The relative phase of local oscillators (LO) and signals (SIG) can be stabilized using 
mirrors mounted on piezo-electric transducers (PZT). The homodyne currents are low-pass hltered and the resulting signals Zd 
and Zr are digitized at 50 MHz. 

(b) Optical micrograph of the mechanical oscillator. A doubly-clamped oscillator from SiN with a Ta 205 /Si 02 Bragg mirror is 
used as the end mirror of the optomechanical cavity. 


b. Amplitude noise 


a. Measurement. We measure the amplitude uoise 
of our laser by direct detectiou ou au luGaAs photodiode 
with a detectiou baudwidth of up to 20 MHz. The statis¬ 
tics of the detected photocurreut I{t) is proportioual to 
the statistics of the uumber operator N = })^h. 

Usiug the defiuitiou (Bl) of the laser amplitude aud 
phase fluctuatious, we hud for the uoise power spectral 
deusity of the photou uumber iu direct detectiou 

Snn{^)= +(/3o'+E[<5/32]), (B2) 

'-^^' 

classical amplitude noise shot noise 


Note that we omitted DC-terms proportioual to 5{(jj) iu 
(B2). Note also that the uoise power spectral deusity of 
the shot uoise is equal to the meau photocurreut. 

Siuce all extra-cavity fields |/3(t)) are derived from the 
same laser field with state 


l7(h) = 


(7o+<57(i))e 


(B3) 


we cau write /3(t) = r^{t) such that /3o = ryo aud 5p{t) = 
r6j{t); here, = P/Piaser is the ratio of the power P of 
the field P{t) we measure to the power Piaser of the laser. 
Therefore, we cau write 


Snn{^) — ^Po^(6/3)(d/3){^) +/^0 

OC 4(P^/Plaser)P((57)((57)(^) + • 

classical amplitude noise shot noise 

For high optical powers P, shot uoise is uegligible aud we 
cau directly measure the classical amplitude uoise, such 
that oc S]sfN{^) ^ By scaliug appro¬ 

priately with the ratio of the powers we cau theu calculate 
the classical amplitude uoise power also for weaker fields. 

b. Modeling. Our goal is to hud a state space model 
that reproduces the spectral depeudeuce of the ampli¬ 
tude uoise i^i ^ satisfactory way. Siuce, for 

high optical powers, the measured photocurreut I{t) is di¬ 
rectly proportioual to the relevaut amplitude uoise sigual 
S/3{t), we use I{t) as iuput for ideutifyiug a suitable state 
space model. This is doue usiug the MATLAB System 
Ideutificatiou toolbox. The results are showu iu Fig. 5. 


c. Homodyne detection with noise 

We use a staudard homodyue setup iu which local os¬ 
cillator (mode 1) aud sigual (mode 2) are combiued ou 
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FIG. 5. Measured and modeled amplitude noise. The mea¬ 
surements (red points) show the photocurrent noise power 
spectral density Sii{(ju) in direct detection for the detuned 
(a) and resonant (b) optical beam with respective shot noise 
(black line) and electronic noise levels (gray line). In the 
relevant frequency range, Sii{(ju) is dominated by the contri¬ 
bution of the classical amplitude noise S(^sp)( 5 p){^)- A state 
space model for the amplitude noise is identified based on 
these measurements. We plot the noise power spectrum of 
an exemplary output time series of these state space models 
(blue points) together with their transfer functions (straight 
lines). 



optomechanical interaction changes the quantum statis¬ 
tics of the output light field, e.g., leading to squeezing 
[45-47]. We then find (in dependence of the homodyning 
angle cp) 

Sn.n. (oj) = {Pl + E[S/3^]){1 + r^) + ^fil] + £[^2] 

shot noise 

+4/^0 {lo) (B5) 

"-V-" 

signal noise 

+(rcos(^)^ 16/?o S'5/36/3(w) 

^ ^ ^ 
classical amplitude noise 

In (B5) we again omitted DC-terms and terms that are of 
more than second order in the fluctuations {6^, (j), /3x, Py). 

Depending on the homodyning angle (^, we detect the 
noise Sx{^)x{(p) (<^) of f he signal in the generalized quadra¬ 
ture x{(p) = cos{(p)^x + sin((/9)/3y. For (/? = 0, in particu¬ 
lar, we detect the amplitude fluctuations Sp^p^{uj) of the 
signal together with the common amplitude fluctuations 
of signal and local oscillator Ssp6p{uj). For (p = 7r/2, 
on the other hand, we detect the phase fluctuations 
S/ 3 yPy{uj). The common phase fluctuations 0(t) of sig¬ 
nal and local oscillator cancel as we aligned the beam 
paths to equal lengths such that both beams have the 
same instantaneous optical frequency at the detectors. 
For any value of the homodyning angle cp, we detect a 
constant shot noise background proportional to the com¬ 
bined optical power of signal and local oscillator field. 


d. Frequency and phase noise 


a 50:50 beam splitter. The output modes 3 and 4 are 
detected and their photocurrents subtracted from each 
other. The statistics of the resulting difference current 
I-{t) is proportional to that of the photon number dif¬ 
ference operator 

N_ = bibs - h\bi = e-^‘^b\b2 + e^n\bi. 

Here, ip is the additional (constant) phase that the local 
oscillator acquires with respect to the signal. 

Now we assume the following states for the local oscil¬ 
lator (mode 1) and signal (mode 2): 

|■*/’2) = |?'(/?0+<^/3)e“*'^+/3a;+«/?y)- (B4) 

The signal is derived from the same laser as the local 
oscillator; therefore, its coherent amplitude has a con¬ 
tribution + 5p)e~'^^ with = Psig/^LO^ the 
ratio between the optical powers Pgig and Ppo of signal 
and local oscillator. The signal beam carries additional 
amplitude and phase quadrature fluctuations j3x and /3y 
due to the interaction with the optomechanical cavity, 
which are assumed to be zero-mean E[Px] = '^[Py] = 0- 
Note that, by writing the signal field state as a displaced 
coherent state (B4), we neglect the possibility that the 


a. Measurement. To measure the phase noise of our 
laser we use delayed self-homodyning [48]. The laser 
beam is split into two beams (beams 1 and 2) of equal 
power. One of the resulting beams (beam 2) is then de¬ 
layed by a time AT before both beams are recombined 
on a 50:50 beam splitter. Both outputs of the 50:50 beam 
splitter are detected and their photocurrents subtracted. 
At low frequencies, where the large-amplitude phase fluc¬ 
tuations occur, the phase between both beams is stabi¬ 
lized to (/2 = 7r/2. 

The difference current /_ is then proportional to the 
number difference operator 

N_ = b\bs - b\bi = i{b\bi - b\b2), 

and we assume beam 1 and 2 to be in the 
states IV^i) = F and |'02) = 

|(/3o + + AT))e“*'^^^+^^^), respectively. 

Assuming small phase fluctuations 1 — i(j){t) 

and neglecting terms of more than second order in the 
fluctuations, we find the noise power spectrum of the 
difference current to be proportional to 

Sn.n. (co) = 4/3^ Sa4>A4>{(^) + 2(/3g + ^[5/3^]). (B6) 
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In (B6), 5 'a 0A(/)(<^) is the noise power spectral density of 
the accumulated phase difference /S.(j)(t) = dr(j){T) 

with 0(r) the instantaneous optical frequency. This can 
be related [49] to the noise power spectrum of the fluc¬ 
tuating phase (p via 

^A0A0(^) 4sin ^ ^ ^ 

such that we get 

Sn.nA'^) = 16/3o sin^ 

classical phase noise contribution 

+2(/32 + E[<5/32]), 

shot noise 

The noise power spectrum of the phase noise 
can therefore be obtained from the spectrum of the de¬ 
tected difference current by subtracting the shotnoise 
and dividing by the function sin^which acts like 
a frequency-dependent gain for the phase fluctuations. 
Note that for frequencies uj around integer multiples of 
27r/AT, the interferometer is insensitive to phase noise. 
To get an accurate phase measurement also at these fre¬ 
quencies, the measurement has to be repeated with a 
different delay AT. 

b. Modeling. We actually need the fluctuating fre¬ 
quency 0 rather than the fluctuating phase 0 as input to 
the optomechanical state space model. To obtain a state 
space model for the frequency noise S^^{uj) from the pho¬ 
tocurrent /_ measured in delayed self-homodyning, we 
note that S. Furthermore, we measure 

at optical powers which are high enough for the shot noise 
to be neglected such that we get 

•2 

Sn{u:)^ -(B7) 

We choose a small AT = 27 ns such that we can linearize 
the sin^-term in (B7) for the relevant frequencies below 
5 MHz and get oc Hence, for high op¬ 

tical powers and low frequencies the measured photocur¬ 
rent /_ is directly proportional to the frequency noise 0 
and can be used as input for identifying a state space 
model in analogy to the amplitude noise (see Fig. 6). 

3. Spurious mechanical modes 

In our experiment, we use a doubly-clamped mechan¬ 
ical oscillator [see Fig. 4(b)], which is coupled to two 
optical fields. This oscillator supports multiple mechani¬ 
cal modes, as is common with clamped oscillators. Some 
mechanical modes of our oscillator are shown in Fig. 7 
which have been calculated by finite element modelling. 

We are interested in coupling to the fundamental out- 
of-plane mechanical mode [Fig. 7(a)], as it exhibits the 



Frequency (Hz) 

FIG. 6. Measured and modeled frequency noise. The 
measurements (red points) show the photocurrent noise 
power spectral density Si_ i_ (a;) obtained using delayed self- 
homodyning for the resonant optical beam. Also shown is 
the electronic noise level (gray line). In the relevant fre¬ 
quency range, is proportional to the classical fre¬ 
quency noise A state space model for the frequency 

noise is identified based on these measurements. We plot the 
noise power spectrum of an exemplary output time series of 
the state space model (blue points) together with its trans¬ 
fer function (straight lines). Note that the sharp resonance 
around 2.5 MHz are resonances of the piezoelectric transducer 
in the laser which we do not model separately. 



FIG. 7. Finite-element modelling of the utilized doubly- 
clamped mechanical oscillator. Shown is the displacement 
amplitude of (a) the fundamental mode with a modelled fre¬ 
quency of 1.28 MHz and higher order modes with frequencies 
of (b) 2.31MHz, (c) 3.08 MHz, (d) 4.31 MHz. 


largest optomechanical coupling strength. However, all 
other mechanical modes will also couple to the optical 
fields, albeit with lower coupling strength. Figure 8 
shows the noise power spectrum of the homodyne signal 
Zd of the detuned beam, where we can attribute some 
peaks in the spectrum to mechanical modes. We also 
incorporate these other mechanical modes in the state 
space model of our experiment. Only then the Kalman 
filter will yield an accurate estimate of the state of the 
fundamental mechanical mode. The consistency of the 
state estimation can be checked by inspecting the inno¬ 
vation sequence, for details see Sec. D. 
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weak coupling (gd=0.2 K) 



FIG. 8. Noise power spectrum of signal Zd in the weak 
coupling regime. The fundamental mechanical mode at a 
frequency of 1.278 MHz is marked with the vertically gray- 
shaded region. Other peaks at frequencies of 2.325 MHz, 
3.05 MHz, 4.237 MHz, 5.604 MHz and 7.24 MHz can be at¬ 
tributed to mechancial modes calculated via finite element 
modelling and are marked at the top horizontal axis, cf. Fig. 7. 


Appendix C: Complete State Space Model 

To construct a state space model for the complete op¬ 
tomechanical experiment, we split it into components of 
independent state space models, which are connected 
in an appropriate way (using network synthesis, see 
[50, 51]). Figure 9 shows a schematic of the complete 
state space model and its components. These compo¬ 
nents describe 

• the evolution and output of the cavity- 
optomechanical system. 


• optical losses. 


• homodyne detection. 


• and various forms of optical noise. 


1. State evolution 

In our experiment, the cavity-optomechanical system 
is comprised of two optical cavity modes driven by two 
external laser fields that interact with a micro-mechanical 
oscillator, which supports multiple mechanical modes. 
We are mainly interested in the coupling of the two intra¬ 
cavity optical fields to the out-of-plane fundamental me¬ 
chanical mode of the mechanical oscillator. Then, the 
state vector consists of the mechanical and optical intra¬ 
cavity quadratures , whose time 

evolution can be described by the linearized Langevin 


equations in the rotating, displaced frame: 

q = UmP (Cla) 

p = -UmQ - ImP + T] ^^(cos 9iXi - sin OiHi) + ^ 

i=r,d 

(Clb) 

Xi = -KXi + Aiyi + gi sindiq + 

+ ■\/2k2x'^2 + + lao.il sin Oi^i 

(Clc) 

ili = -KUi - AiXi + Qi cos Oiq + 

+ y/2K2y'^2 + I *^0,11 COS Oi^i 

(Cld) 

The time evolution of the state vector is driven by 
(white) Brownian thermal noise from the mechanical 
bath as well as by noise on the driving laser fields, which 
consists of classical amplitude noise dp{t) and phase noise 
0(t), as well as shot noise, which is introduced via the 
terms with i G {d^r} for the detuned and reso¬ 

nant beam and j G {1,2} for the input coupling mirror 
and end mirror. 

Shot noise and mechanical thermal noise are white and 
are incorporated in the state space model directly as 
process noise and, in the case of measurement shot 
noise, as measurement noise v^. However, classical laser 
noise in our experiment is not white. We model it inde¬ 
pendently as described in Sec. A 2 by constructing state 
space models, whose output resembles the measured ex¬ 
perimental noise characteristics. The output vectors of 
the noise state space models are treated as deterministic 
input Ut = {5pd{i)i d/3r{t), 0^(t))^, which enter the 

dynamic and measurement equation of the optomechan¬ 
ical system. 

The matrices governing the evolution of the state vec¬ 
tor are the process matrix A G the input matrix 

B G the noise matrix L G and the process 

noise covariance matrix W G The process matrix 

A describes the evolution of the optomechanical system 
and can be directly read off from the Langevin equations: 

(A^i^gd-jOd) 

A= A 4 (gd,dd) A 2 (Ad) 0 . (C2) 

\ 0 A 2 {Ar) J 

Note that matrix (C2) and the following matrices are 
written in block form. Thus, zeros stand for zero matrices 
with appropriate dimensions. In (C2), the evolution of 
the mechanical and optical quadratures is described by 
the matrices 

Miu,!) = , (C3) 

M^) = (Zl _t) , (C4) 
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where uJm and 7^ are the frequency and linewidth 
(FWHM) of the mechanical mode and A and k the detun¬ 
ing and linewidth (HWHM) of the optical cavity mode, 
respectively. The interaction between mechanical and 
optical modes, on the other hand, is described by the 
matrices 

A 3 {g,d) =( ^ n d (C5) 

^ ygcosO —gsmOJ ^ ’ 

(“) 

where g = y/2agQ is the linearized optomechanical 
coupling strength and the angle 0 = arctan(A//^) 

parametrizes the detuning from cavity resonance. 

To accurately describe the experimental situation we 
actually have to model multiple mechanical modes of the 
mechanical oscillator. The generalization of matrix (C2) 
to the multimode case is straightforward: for each addi¬ 
tional mechanical mode /c, one determines the appropri- 
ate matrices 

from the corresponding optomechanical parameters and 
adds them as new blocks to (C2). Note that adding M 
additional mechanical modes increases the state space di¬ 
mension from 6 X 6 to (6 + 2M) x (6 + 2M). 

The input matrix B determines the coupling of the 
classical amplitude and phase noise of the laser to the 
optical intra-cavity quadratures 


/ 0 0 \ 

Bi(0d,ao,d) 0 , (C7) 

y 0 Bi{0r,aQ^r)/ 

with 


Bi{0, 0^0,i) 


/x/2^ |(ao,i| sin6>\ 
y 0 |<ao,i| cosOJ ’ 


(C8) 


where k,i is the optical decay rate through the input cou¬ 
pler mirror and |ao,i| = \J J 1)^2 yj the intra-cavity 

photon number. These matrices can again be directly 
inferred from the Langevin equations, when optical noise 
is included in their derivation (see, e. g., [34]). 

White thermal noise and shotnoise, on the other hand, 
are treated as process noise w^. Its coupling to the state 
vector evolution is given by the noise matrix 


(Li{^m) 0 0 0 0 \ 

L=\ 0 L2{i^i) 12 ( 1 ^ 2 ) 0 0 , (C9) 

V 0 0 0 L2{ki) L2{k2)J 

with 

Li(7) = -V^(^jy (CIO) 


Li{j) describes the driving of the mechanical mode by 
white thermal noise (assumed to act only on the mechan¬ 
ical momentum quadrature) and 


L2{i^) = -V^i, 


describes driving of the optical modes by shot noise (with 
1 as identity matrix of appropriate dimensions). Note 
that in our case, shot noise enters both from mirror 1 
(input coupler) and mirror 2 (mechanical oscillator) with 
rates ki and K 2 = respectively, where K 2 also 

incorporates additional intra-cavity loss. This accounts 
for the double-sidedness of the cavity. 

The process noise covariance matrix is 

W = diag(n^ + 1/2, V2 ,..., V2), (C12) 

where rim = i^m{T^uJm) ~ ki)T/{huJm) denotes the ther¬ 
mal occupation of the mechanical bath and 1/2 is the 
shot noise contribution of the driving laser fields. 


2. Measurement 

To correctly describe the effect of noise and losses on 
the measurement, we formally split the measurement into 
several steps: 

1. the intra-cavity quadratures {xd^Vd^Xr^Vr) are re¬ 

lated to the quadratures of the output fields 
{x ^^^, ) via cavity input-output rela¬ 

tions; the output quadratures are affected by clas¬ 
sical and quantum laser noise, 

2. a state-space model for the optical loss; this 

takes the output quadratures ^/^^) 

as input and returns the attenuated quadratures 

that arrive at the homodyne 

detectors, 

3. a measurement of one (generalized) quadrature per 
returning laser field, which yields the actual mea¬ 
surements {Zd^Zr). 

In conjunction, these steps guarantee that the measure¬ 
ments are physical (i. e. respect uncertainty relations) 
and that all noise sources are correctly taken into ac¬ 
count. Note, we also model the detector spectral response 
as a digital filter through a separate state space model 
(see Sec. A 2). 


a. Cavity output 


The matrices relating the quadratures of the intra¬ 
cavity laser fields to those of the extra-cavity laser fields 
are the measurement matrix C G the through¬ 
put matrix D G the measurement noise covariance 

matrix V G and the noise cross-correlation matrix 

M G 

The measurement matrix is 


with 


fO CM 0 \ 

1^0 0 


(C13) 


(Cll) 


Ci(k) = v^l. 


(C14) 
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In (C13), ni is the decay rate of the input coupler mirror 
because we only measure light coming from the input 
coupler. 

Further, the throughput matrix is 

® = (^o' D.) 

with 

D. = (J J) . (C16) 

Remember that D in (Alb) describes the feedthrough of 
a deterministic input to the state-space model onto the 
measurement, which is in our case classical laser noise. 
Hence, D describes the impact of incoming classical laser 
noise on the extra-cavity quadratures of the returning 
laser fields. 

The fact that (C16) has an entry only at position (1,1) 
means that effectively only amplitude noise is reflected off 
the cavity. In principle, also phase noise ^i{t) is reflected 
off the cavity; but we can neglect it here since it cancels 
in homodyne detection (assuming that the local oscillator 
and reflected optical field have the same instantaneous 
optical frequency). Note that phase noise does, however, 
affect the intra-cavity optomechanical state evolution (as 
described previously). 

The measurement noise covariance matrix is 


with 

Dc{t) = cosrU, (C21) 

T),(t) = sinrU (C22) 

and r = arccos (y^). In our case, (xi,^i) are the extra¬ 
cavity quadratures and (^2,^2) are quadra¬ 

tures of optical shotnoise. The output quadratures are 
denoted as 

The optical shotnoise is generated by a separate state 
space model, which is only made up of the measurement 
noise covariance matrix 

1 / = il. (C23) 

Note that each of the two cavity-output beams is fed into 
a separate loss model with individual loss parameters rjd 
and r]r. 

c. Homodyne detection 

The state space model for homodyne detection guar¬ 
antees the physicality of the measurement. It maps 
to the actual physical measure¬ 
ments Zt using the following throughput matrix Dt G 
M2x4. 


R = diag(V2,V2,V2,V2), (CI7) 


i.e., shot noise. Finally, the cross-correlation matrix be¬ 
tween process and measurement noise is 


M = 


Ml 


0 

V 0 



(C18) 


with 

Ml = b. (C19) 

The matrix M guarantees that the optical shot noise con¬ 
tribution to the process and measurement noise is per¬ 
fectly correlated. 


with 


A 


0 DiiiPrit))J ^ 


Di{ip) = (cosi^ sini^) . 


(C24) 


(C25) 


The measurement zj = {zd{ipd{t)),Zr{‘fir{t)))'^ can be 
(explicitly) time-dependent if the phase ^pi (t) of the local 
oscillator changes in time. 

Note that we do not have to include noise in this final 
measurement step since all optical noise (shot noise as 
well as classical noise) has already been included in the 
quadratures ). 


Appendix D: Kalman filter consistency 


b. Optical loss 


1. Statistical consistency tests 


Optical loss is commonly modeled as a beam splitter 
with intensity transmission r] [52], where the signal en¬ 
ters through one port (port I) and optical shot noise 
through the other (port 2). The state space model of 
the beam splitter is then only given by the through¬ 
put matrix D G which maps the input quadra¬ 

tures u = (^1,^1,^2,^2)^ to the output quadratures 
Z = (3:3,2/3,3:4,2/4)'^: 


D = 


( D^t) D,{t)\ 


(C20) 


The Kalman filter must be checked for a correct oper¬ 
ation in two ways [41]. First, we have to test the filter 
consistency for a given model, i.e., that the estimate x 
converges to the true state x. In the case of the Kalman 
filter this means that the sample mean of the estima¬ 
tion error = x^ — Xt goes to zero (which means that 
the filter is unbiased) and the respective sample covari¬ 
ance matrix converges to Pf Typically such tests can 
only be implemented in simulations, as the true state 
is normally not available in an experiment. Second, we 
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FIG. 9. Schematic overview of the complete system model for describing the cavity-optomechanical experiment. State space 
models of optical noise (Lorentz noise, colored noise), the cavity-optomechanical system, optical loss and homodyne detection 
are combined via network synthesis to form a single state space model for modeling the total experiment. 


have to ensure that the state-space model accurately de¬ 
scribes the experimental system, which means that all 
experimental parameters have been determined correctly 
and possible approximations of the system dynamics (in¬ 
cluding the description of all noise processes) are fulfilled 
to a sufficient degree. The so-called innovation sequence 
Ut = Zt — Ctkt —which records the difference between 
the predicted and the actual measurement results—is our 
main tool to test the consistency of the Kalman filter 
with the real experimental situation. For a model per¬ 
fectly matching the experimental system the innovation 
sequence is a zero-mean white Gaussian process with a 
covariance matrix given by St = CtPtCj S V [41, 42]. 
For convenience we introduce the normalized process 

Dt = (Dla) 

(where St = LtLj is the Cholesky decomposition of St) 
whose components now have unit variance. Both the 
whiteness of the process and the correct distribution can 
be tested for an experimentally recorded innovation se¬ 
quence. In order to test for whiteness we calculate the 
so-called periodogram^ of a recorded sequence, which 
asymptotically must be distributed with two degrees 
of freedom. The Gaussianity can be tested by calculating 
the sample’s cummulative distribution function (CDF) or 
the corresponding probability density function (PDF). 


^ The periodogram at Fourier frequency fj = j/N is defined as 
Pfj) — with the Fourier coefficients Xc{fj) = 

l/\^EtliCOs(27r/jt)a:t, Xs{fj) = sin (27r/jt)a:t, 

N total samples and the random process x with values xt at 
time t. Hence, the periodogram is a sum of independent, squared 
zero-mean, white Gaussian random vectors and as such asymp¬ 
totically X^(2) distributed with 2 degrees of freedom. 


weak coupling strong coupling 

sample mean of i^d 

0.004 

-0.012 

sample mean of i^r 

0.031 

0.031 


TABLE 1. Mean values of normalized innovations for the weak 
and strong coupling regime. The expected value is zero. 


2. Consistency test results 

Here we show the results of the statistical consistency 
analysis for the data presented in the main text. Table 

1 lists the mean values of the components Ud and rep¬ 
resenting the innovation sequence for the measurements 
of the detuned and resonant beam, respectively. We find 
that the mean of the normalized innovations is close to 
zero. To test the distribution of i>d and i>r we calcu¬ 
late their conditional distribution functions and proba¬ 
bility density functions, see Fig. 10(a). We find that they 
closely match the expected Gaussian distribution. To be 
quantitative, we compute the ratio of i>i (i G d,r) that 
are expected to lie within a two-sided 95% confidence re¬ 
gion (i. e., between 0.025% and 0.975%) of the expected 
Gaussian distribution. For our data, we find values close 
to the expected 95%. 

The small deviation can be more closely inspected by 
considering the noise power spectrum of Pt, which is 
shown in Fig. 10(b). We smooth the periodogram us¬ 
ing Welch’s method using a Hamming window and split¬ 
ting the data into eight non-overlapping segments. Then 
the resulting noise power spectrum is distributed with 

2 • 8 = 16 degrees of freedom. In Fig. 10(b) we plot 
the two-sided 95% confidence region of this distribution 
as horizontal gray region behind the data. We see that 
most datapoints lie within this region and thus can be 
considered white. Most deviations are observed in the 
low frequency regime around 200 kHz (a dip in the noise 
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(a) weak coupling (gd=0.2 k) 



strong coupling (gd=1.68 k) 
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FIG. 10. Statistical analysis of normalized innovation se¬ 
quence. The cumulative distribution function (CDF) and 
probability density function (PDF) of the normalized innova¬ 
tions Vi are shown in (a) for the detuned (red) and resonant 
(blue) beam, each for the weak (left) and strong (right) cou¬ 
pling regime. The normalized innovations are expected to be 
a zero-mean Gaussian white noise process with unit variance, 
i.e., Vi ~ A/'(0,1). The numbers in the CDF panel represent 
the fraction of Vi that are contained in the expected two-sided 
95% confidence region (±2cr) of a random variable distributed 
as A/’(0,1). The noise power spectrum, Snmi, of Vi is shown 
in (b). The horizontal, gray shaded area indicates the 95% 
confidence region of the expected x^(16) distribution. The 
numbers in the panels reflect the actual value of experimental 
Vi contained in this region. The vertical, gray shaded area in¬ 
dicates the region, where the mechanical mode contains most 
of its noise power. 


The resolution of the digital calculations needs to be 
much greater than the minimal step size of the harmonic 
evolution of an oscillation, i.e., the amplitude of a single 
phonon. This enables the possible operation at the stan¬ 
dard quantum limit. Second, the dynamic range must 
cover the thermal state of the oscillator. Assuming a 
typical FPGA clock rate of 500 MHz, and the parame¬ 
ters of the mechanical mode of interest, we obtain a nec¬ 
essary resolution of fmdt/y/rith ~ 4 • 10' -7 ^ 2-2% This 
dynamic range can be covered by a 27 bit fix-point cal¬ 
culation, as used in the StratixV GS FPGA. This FPGA 
chip can process 1590 multiplications per clock cycle. It 
can be used, e. g., to calculate a steady-state state space 
model of one mechanical mode at 1.27 MHz that is cou¬ 
pled to two optical modes, and incorporates 25 additional 
mechanical modes and 10 dimensions for modeling col¬ 
ored laser noise. A final criterion concerns the delay due 
to the detection and the processing of the signals, which 
needs to be taken into account. We obtain about 50 ns 
system delay, when accounting for 4 m of optical path 
length and cables, the delay from signal converters and 
the processing in a 500 MHz FPGA. This amounts to 
a phase shift of 23° at 1.27 MHz, the frequency of the 
mechanical mode of interest. Therefore efficient feed¬ 
back control of mechanical oscillators with a resonance 
frequency on the order of 1 MHz is feasible with state-of- 
the-art technology. 


power) and in the frequency band 2 — 5 MHz (many sharp 
peaks). Both features are not considered in the model for 
the Kalman filter and therefore lie outside the 95% confi¬ 
dence region. The series of sharp peaks can be attributed 
to resonances of a piezoelectric transducer, which is at¬ 
tached to the laser crystal of the driving laser. The spec¬ 
tral dip at 200 kHz we attribute to the spectral response 
of the photo detectors. 

Finally we note that we also tested the consistency 
of the implemented Kalman filter for simulations of a 
perfectly matched system and find excellent agreement 
with the criteria presented above. 


Appendix E: Real-time feedback 

The execution of the Kalman filter requires calcula¬ 
tions of extensive matrix multiplications with a specific 
dynamic range and minimal delay. Therefore, the imple¬ 
mentation of a real-time Kalman filter hinges on a fast 
and precise hardware platform, which can be based on a 
high resolution, field programmable gate array (FPGA). 
Such a hardware platform has to fulfill several criteria. 

First, a minimum precision level has to be achieved. 
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Appendix F: Symbols and definitions 


symbol 

meaning 

xt, x(t) 

state vector (random process) 

Zt 

measurement vector (random process) 

Wt 

process noise 

Vt 

measurement noise 

Xt, x{t) 

estimate for xt 

Sxx (^) 

noise power spectrum of Xt 

(•) 

quantum expectation value 

E[.] 

classical (population) expectation value 

At 

process matrix 

Bt 

input matrix 

Ct 

measurement matrix 

Dt 

throughput matrix 

W 

process noise covariance matrix 

V 

measurement noise covariance matrix 

M 

noise cross-correlation matrix 

Kt 

Kalman gain matrix 

Ft 

estimation error covariance matrix 


TABLE IL State space model specific symbols. 


symbol definition 

name 

value 

LOm 


mechanical frequency 

27r • 1.278 • 10® Hz 



mechanical linewidth (full width at half maximum, EWHM) 

27r • 265 Hz 

K 

Kl + 1^2 

cavity decay rate (half width at half maximum, HWHM) 

0.341wm 

1^1 


input-coupler cavity decay rate (HWHM) 

0.2775(x;m 

1^2 


total output-coupler and loss cavity decay rate (HWHM) 

0.0635cJrn 

9d 

V^\ao^d\9o 

optomechanical coupling rate detuned beam 

between 0.2/^ and 1.68/^ 

Qr 

A/2|Q^0,r 

optomechanical coupling rate resonant beam 

0.2/^ 

90,i 

<^o,i / h 

single-photon coupling rate 

27r-7.7Hz 

L Y rneffo;^ 

arctan (^) 

Oi 



Ai 

effective detuning 

At- — 0, A(^ — OJrti 

Ao,z 

Uc — 

cavity-laser detuning 


OJc 


cavity resonance frequency 



2i^cl\i 

laser frequency 

~ 1.77- 10^® Hz 

Pi 


extra-cavity optical power 


1^0,i 

Y ^^0 

extra-cavity amplitude 




classical amplitude fluctuation of laser drive (extra-cavity) 




classical phase noise of laser drive (intra-cavity and extra-cavity) 




intra-cavity steady-state amplitude 


Sai{t) 


classical amplitude fluctuation of laser drive (intra-cavity) 


(Pi 


phase between local oscillator and signal beam in homodyne detection 



TABLE III. Optomechanical symbols and definitions. The subscripts i G {d, r} denote the two laser beams. 



16 


[1] Robert F Stengel, Optimal Control and Estimation 
(Dover Publications, Mineola, NY, 1994). 

[2] H Wiseman and Gerard J Milburn, Quantum Measure¬ 
ment and Control (Cambridge University Press, Cam¬ 
bridge, 2010). 

[3] RE Kalman, “A New Approach to Linear Filtering and 
Prediction Problems,” Journal of Fluids Engineering , 
35-45 (1960). 

[4] R E Kalman and R S Bucy, “New Results in Linear Fil¬ 
tering and Prediction Theory,” Journal of Fluids Engi¬ 
neering 83, 95-108 (1961). 

[5] MS Grewal and AP Andrews, “Applications of Kalman 
Filtering in Aerospace 1960 to the Present [Historical 
Perspectives],” IEEE Control Systems 30, 69-78 (2010). 

[6] Lee Samuel Finn and Soma Mukherjee, “Data condition¬ 
ing for gravitational wave detectors: A Kalman filter for 
regressing suspension violin modes,” Physical Review D 
63, 1-19 (2001). 

[7] JM Geremia, John K Stockton, Andrew C Doherty, and 
Hideo Mabuchi, “Quantum Kalman Filtering and the 
Heisenberg Limit in Atomic Magnetometry,” Physical 
Review Letters 91, 250801 (2003). 

[8] Hidehiro Yonezawa, Daisuke Nakane, Trevor A Wheat- 
ley, Kohjiro Iwasawa, Shuntaro Takeda, Hajime Arao, 
Kentaro Ohki, Koji Tsumura, Dominic W Berry, Timo¬ 
thy C Ralph, Howard M Wiseman, Elanor H Hunting- 
ton, and Akira Furusawa, “Quantum-Enhanced Optical- 
Phase Tracking,” Science 337, 1514-1517 (2012). 

[9] Mankei Tsang, “Time-symmetric quantum theory of 
smoothing,” Physical Review Letters 102 , 250403 (2009). 

[10] Mankei Tsang, Jeffrey H Shapiro, and Seth Lloyd, 
“Quantum theory of optical temporal phase and instan¬ 
taneous frequency. H. Continuous-time limit and state- 
variable approach to phase-locked loop design,” Physical 
Review A - Atomic, Molecular, and Optical Physics 79, 
053843 (2009). 

[11] T A Wheatley, D W Berry, H Yonezawa, D Nakane, 
H Arao, D T Pope, T C Ralph, H M Wiseman, A Fu¬ 
rusawa, and E H Huntington, “Adaptive optical phase 
estimation using time-symmetric quantum smoothing,” 
Physical Review Letters 104, 5-8 (2010). 

[12] V B Braginsky and A B Manukin, “Ponderomotive effects 
of electromagnetic radiation,” Soviet Physics JETP 25, 
653-655 (1967). 

[13] Vladimir B Braginsky and Farid Ya Khalili, Quantum 
Measurement (Cambridge University Press, 1995). 

[14] Markus Aspelmeyer, Tobias J Kippenberg, and Florian 
Marquardt, “Cavity optomechanics,” Reviews of Modern 
Physics 86, 1391 (2014). 

[15] M Paternostro, S Gigan, M S Kim, F Blaser, H R Bohm, 
and M Aspelmeyer, “Reconstructing the dynamics of a 
movable mirror in a detuned optical cavity,” New Journal 
of Physics 8, 107 (2006). 

[16] T A Palomaki, J W Harlow, J D Teufel, R W Simmonds, 
and K W Lehnert, “Coherent state transfer between itin¬ 
erant microwave fields and a mechanical oscillator.” Na¬ 
ture 495, 210-4 (2013). 

[17] D Rugar and P Griitter, “Mechanical parametric amplifi¬ 
cation and thermomechanical noise squeezing,” Physical 
Review Letters 67, 699-702 (1991). 

[18] Y Hadjar, P F Cohadon, C G Aminoff, M Pinard, 


and A Heidmann, “High-sensitivity optical measurement 
of mechanical Brownian motion,” Europhysics Letters 
(EPL) 47, 545-551 (1999). 

[19] T Briant, P-F Cohadon, M Pinard, and A Heidmann, 
“Optical phase-space reconstruction of mirror motion at 
the attometer level,” The European Physical Journal D 
22 , 131-140 (2003). 

[20] Kohjiro Iwasawa, Kenzo Makino, Hidehiro Yonezawa, 
Mankei Tsang, Aleksandar Davidovic, Elanor Hunting- 
ton, and Akira Furusawa, “Quantum-Limited Mirror- 
Motion Estimation,” Physical Review Letters 111, 
163602 (2013). 

[21] D J Wilson, V Sudhir, N Piro, R Schilling, A Ghadimi, 
and T J Kippenberg, “Measurement and control of a 
mechanical oscillator at its thermal decoherence rate,” 
(2014), arXiv:1410.6191v2. 

[22] Howard Carmichael, An open systems approaeh to quan¬ 
tum opties (Springer-Verlag, 1993). 

[23] Crispin W Gardiner and Peter Zoller, Quantum noise, 
3rd ed. (Springer, 2004). 

[24] Luc Bouten, Ramon Van Handel, and Matthew R James, 
“An Introduction to Quantum Filtering,” SIAM Journal 
on Control and Optimization 46, 2199-2241 (2007). 

[25] Stefano Mancini, David Vitali, and Paolo Tombesi, “Op¬ 
tomechanical Cooling of a Macroscopic Oscillator by Ho¬ 
modyne Feedback,” Physical Review Letters 80, 688-691 
(1998). 

[26] A C Doherty and K Jacobs, “Feedback control of quan¬ 
tum systems using continuous state estimation,” Physical 
Review A 60, 2700-2711 (1999). 

[27] Asa Hopkins, Kurt Jacobs, Salman Habib, and Keith 
Schwab, “Feedback cooling of a nanomechanical res¬ 
onator,” Physical Review B 68, 235328 (2003). 

[28] Stefan Danilishin, Helge Muller-Ebhardt, Henning Re- 
hbein, Kentaro Somiya, Roman Schnabel, Karsten Danz- 
mann, Thomas Corbitt, Christopher Wipf, Nergis Maval- 
vala, and Yanbei Chen, “Creation of a quantum oscillator 
by classical control,” 0809.2024 (2008). 

[29] Helge Miiller-Ebhardt, Henning Rehbein, Chao Li, Ya- 
sushi Mino, Kentaro Somiya, Roman Schnabel, Karsten 
Danzmann, and Yanbei Chen, “Quantum-state prepara¬ 
tion and macroscopic entanglement in gravitational-wave 
detectors,” Physical Review A 80, 043802 (2009). 

[30] VP Belavkin, “Optimal filtering of Markov signals with 
quantum white noise,” Radio Eng Electron Physics 25, 
1445-1453 (1980). 

[31] Sebastian G. Hofer and Klemens Hammerer, 

“Entanglement-enhanced time-continuous quantum 

control in optomechanics,” Physical Review A 91, 
033822 (2015). 

[32] Vittorio Giovannetti and David Vitali, “Phase-noise mea¬ 
surement in a cavity with a movable mirror undergo¬ 
ing quantum Brownian motion,” Physical Review A 63, 
023812 (2001). 

[33] P Rabl, C Genes, K Hammerer, and M Aspelmeyer, 
“Phase-noise induced limitations on cooling and coherent 
evolution in optomechanical systems,” Physical Review A 
80, 063819 (2009). 

[34] M Abdi, Sh Barzanjeh, P Tombesi, and D Vitali, “Effect 
of phase noise on the generation of stationary entangle¬ 
ment in cavity optomechanics,” Physical Review A 84, 



17 


032325 (2011). 

[35] R Ghobadi, A R Bahrampour, and C Simon, “Op¬ 
tomechanical entanglement in the presence of laser phase 
noise,” Physical Review A 84, 063827 (2011). 

[36] C Genes, D Vitali, P Tombesi, S Gigan, and M As- 
pelmeyer, “Ground-state cooling of a micromechanical 
oscillator: Gomparing cold damping and cavity-assisted 
cooling schemes,” Physical Review A 77, 033804-9 
(2008). 

[37] Ferdinand Brennecke, Stephan Ritter, Tobias Donner, 
and Tilman Esslinger, “Gavity optomechanics with a 
bose-einstein condensate,” Science 322, 235-238 (2008). 

[38] Kater W March, Kevin L Moore, Subhadeep Gupta, 
and Dan M Stamper-Kurn, “Observation of quantum- 
measurement backaction with an ultracold atomic gas,” 
Nature Physics 4, 561-564 (2008). 

[39] T P Purdy, R W Peterson, and G A Regal, “Observa¬ 
tion of Radiation Pressure Shot Noise on a Macroscopic 
Object,” Science 339, 801-804 (2013). 

[40] Helge Miiller-Ebhardt, Henning Rehbein, Roman Schn¬ 
abel, Karsten Danzmann, and Yanbei Ghen, “Entan¬ 
glement of Macroscopic Test Masses and the Standard 
Quantum Limit in Laser Interferometry,” Physical Re¬ 
view Letters 100, 013601 (2008). 

[41] Yaakov Bar-Shalom, X Rong Li, and Thiagalingam 
Kirubarajan, Estimation with Applications to Tracking 
and Navigation: Theory Algorithms and Software (John 
Wiley & Sons, 2001). 

[42] Eerdinand van der Heijden, Robert Duin, Dick de Rid- 
der, and David M J Tax, Classification, Parameter Esti¬ 
mation and State Estimation: An Engineering Approach 
Using MATLAB (John Wiley & Sons, 2005). 

[43] Robert E Stengel, Optimal Control and Estimation, reis¬ 
sue edition ed. (Dover Publications, New York, 1994). 

[44] A A Glerk, M H Devoret, S M Girvin, Elorian Mar- 
quardt, and R J Schoelkopf, “Introduction to quantum 


noise, measurement, and amplification,” Reviews of Mod¬ 
ern Physics 82, 1155-1208 (2010). 

[45] Daniel W G Brooks, Thierry Botter, Sydney Schreppler, 
Thomas P Purdy, Nathan Brahms, and Dan M Stamper- 
Kurn, “Non-classical light generated by quantum-noise- 
driven cavity optomechanics,” Nature 488, 476-480 
( 2012 ). 

[46] Amir H Safavi-Naeini, Simon Groblacher, Jeff T Hill, 
Jasper Ghan, Markus Aspelmeyer, and Oskar Painter, 
“Squeezed light from a silicon micromechanical res¬ 
onator,” Nature 500, 185-189 (2013). 

[47] T P Purdy, P-L Yu, R W Peterson, N S Kampel, and 
G A Regal, “Strong optomechanical squeezing of light,” 
Physical Review X 3, 031012 (2013). 

[48] Amir H Safavi-Naeini, Jasper Ghan, Jeff T Hill, Simon 
Groblacher, Haixing Miao, Yanbei Ghen, Markus As¬ 
pelmeyer, and Oskar Painter, “Laser noise in cavity- 
optomechanical cooling and thermometry,” New Journal 
of Physics 15, 035007 (2013). 

[49] Ralf Riedinger, Optomechanical State Reconstruction and 
Optical Noise Reduetion for Cavity Optomechanies Ex¬ 
periments, Master thesis, Philipps University Marburg, 
Germany (2013). 

[50] J Gough and MR James, “The Series Product and Its 
Application to Quantum Eeedforward and Eeedback Net¬ 
works,” Automatic Gontrol, IEEE Transactions on 54, 
2530 -2544 (2009). 

[51] H Nurdin, M James, and A Doherty, “Network Synthe¬ 
sis of Linear Dynamical Quantum Stochastic Systems,” 
SIAM Journal on Gontrol and Optimization 48, 2686- 
2718 (2009). 

[52] Ulf Leonhardt, “Quantum statistics of a lossless beam 
splitter: SU(2) symmetry in phase space,” Physical Re¬ 
view A 48, 3265-3277 (1993). 


